Vortex liquid crystals in anisotropic type II superconductors 
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In a type II superconductor in a moderate magnetic field, the superconductor to normal state 
transition may be described as a phase transition in which the vortex lattice melts into a liquid. In 
a biaxial superconductor, or even a uniaxial superconductor with magnetic field oriented perpendic- 
ular to the symmetry axis, the vortices acquire elongated cross sections and interactions. Systems 
of anisotropic, interacting constituents generally exhibit liquid crystalline phases. We examine the 
possibility of a two step melting in homogeneous type II superconductors with anisotropic super- 
fluid stiffness from a vortex lattice into first a vortex smectic and then a vortex nematic at high 
temperature and magnetic field. We find that fluctuations of the ordered phase favor an instability 
to an intermediate smectic-A in the absence of intrinsic pinning. 



Recently, there has been much interest generated con- 
cerning high temperature superconductors in a magnetic 
field. Various experiments have studied the interplay of 
superconductivity with coexistent magnetic orders in the 
presence of an external field. Experiments using both 
neutron scattering[|[ || and STM||^ show that there is 
significant local electronic inhomogeneity. These and 
other experiments lend credence to the idea that there 
may be electronic liquid crystalline phases in strongly 
correlated systems, leading to anisotropy even within a 
Cu02 plane. 

The cuprate superconductors are also ideal laborato- 
ries for studying vortex physics, due to the large values 
of K = Xab/iab (whcrc \ab and S^ab are the London pen- 
etration depth within a plane, and the coherence length 
within a plane, respectively) and small critical depinning 
current. Q In this letter, we consider the effects of an 
anisotropic superfluid stiffness on the vortex phases in 
superconductors in the continuum limit. We account for 
this anisotropy by allowing different effective masses in 
the three crystalline directions, which we will call TOq, 
TTifc, and TOc. 

In a biaxial superconductor (with different effective 
mass in each crystalline direction), or a uniaxial su- 
perconductor (the effective mass differs in one direction 
only) with magnetic field oriented perpendicular to the 
symmetry axis, vortices acquire elliptical cross sections, 
whether measured by the shape of the core, or by the 
profile of the screening currents or magnetic field which 
penetrates beyond the core. Systems of anisotropic inter- 
acting constituents generically lead to liquid crystalline 
phases. In such a system, we expect the melting to pro- 
ceed from the body centered rectangular lattice, to a 
smectic, to a nematic, as in Fig. |l| (In this case, the 
high temperature phase is trivially nematic due to the 
explicitly broken rotational symmetry introduced by the 
mass anisotropy.) 

Liquid crystals lie somewhere between the full trans- 
lational and rotational symmetry of a liquid, and that 
of a 3D crystal, which has broken rotational symmetry. 
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FIG. 1: Schematic phase diagram of vortex matter with 
anisotropic interactions. Solid lines represent phase transi- 
tions, and the dotted line represents the crossover at Hc2. 
There may also be melted or partially melted phases near 
Hci, between the regions marked Meissner and Lattice. 



and retains only discrete translational symmetry in the 
three directions of the crystal axes. In a superconduc- 
tor, the application of an external magnetic field to pro- 
duce vortices explicitly breaks rotational symmetry. We 
choose axes such that BW'z. In the vortex system, smectic 
phases correspond to liquid-like correlations (and unbro- 
ken translational symmetry) in one direction in the xy 
plane, and simultaneous solid-like correlations (and only 
discrete translational symetry) in the other direction in 
the xy plane. Smectics may be further classified by which 
direction the "elongated molecule" is pointing on average 
with respect to the orientation of the liquid-like layers. 
Call 9 the angle that the long axis of the molecule makes 
with respect to the normal of the liquid- like layers. For 
9 ^ {], the phase is smectic-A, illustrated in Fig. ^ For 
all other values of 9 the phase is a smectic-C. A ne- 
matic phase is characterized by unbroken translational 
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FIG. 2: Smectic-A. Cross sections of vortices in an anisotropic 
superconductor are represented schematically above as filled 
ellipses, for a magnetic field oriented perpendicular to the 
page. There may be an intermediate melting from the elon- 
gated lattice first to a smectic-yl, shown above, and then to a 
nematic phase at high temperature. 



symmetry, with broken orientational symmetry. Corre- 
lations in this phase are liquid-like in all directions, but 
the constituent molecules have a preferred orientation. 

When the magnetic field is oriented parallel to the 
planes 1^, ^, Q in a layered superconductor, the explicit 
translational symmetry breaking of the planes may cause 
the vortex lattice to melt first along the direction of 
the planes, leading to a smectic-C with 9 = ^. Smec- 
tic phases have also been predicted in the presence of a 
driving current [^, as well as chain states which may arise 
when the field is tilted sufficiently away from a crystalline 
Here we explore the possibility of liquid crys- 
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talline phases due to explicit rotational symmetry break- 
ing (mass anisotropy), with no explicit translational sym- 
metry breaking in the problem {i.e. no intrinsic pinning). 
We find that thermal fluctuations of such an anisotropic 
vortex lattice favor an instability to a vortex smectic-A. 
Such a phase, if it exists, should be detectable by several 
experimental probes, including Bitter decoration experi- 
ments, neutron scattering, /xSR, and resistivity measure- 
ments. 

We use continuum elasticity theory to describe the 
thermal fluctuations of an ordered vortex lattice in three 
dimensions. We choose axes such that the external 
magnetic field B\\z. The displacement of a vortex from 
its equilibrium position is denoted by the vector u = 
{ux,Uy) , which is a function of the position z along a 
given vortex. To second order in u, the free energy is 
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where the matrix C contains the elastic constants:! 
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Here we orient the magnetic field perpendicular to the 
axis of symmetry: ~§ _L c| |z. We have also assumed a uni- 
axial superconductor, — mi, = niab Wc ; the elastic 
constants are not yet known for the fully anisotropic, bi- 
axial case. Nonetheless, this uniaxial geometry captures 
the physics we are interested in, namely anisotropic inter- 
actions. We use elastic constants derived from Ginzburg- 
Landau theory. |n|, |l^ The bulk modulus, cii(k), de- 
scribes the compressibility of the lattice. The hard tilt 
modulus, C44(k), corresponds to tilts along the symmetry 
axis c, and the (smaller) easy tilt modulus, c|4(k), corre- 
sponds to tilts perpendicular to c. Similarly, it is easier to 
shear vortices (c|g) along the major axis of the cross sec- 
tional ellipse, rather than perpendicular to it (cgg). Note 
that since the magnetic field is oriented along a crystal 
symmetry axis, there is no mixing between the bulk, tilt, 
and shear moduli. The bulk and tilt moduli are highly 
momentum dependent. In fact, the two soften signifi- 
cantly at the Brillouin zone edge, so that their momen- 
tum dependence is important in the physics of melting. 
The shear moduli are approximately independent of the 
wavelength of the distortion, and we neglect their weak 
momentum dependence. 

We use an extension of the Lindcmann criterion to the 
case of anisotropy]^ and allow for the possibility that 
the lattice may melt in one direction before the other. 
In this case, fluctuations in the x direction compete with 
the lattice spacing in the x direction, and fluctuations in 
the y direction with the lattice spacing in the y direction: 
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is the lattice spacing for the 



triangular lattice of the isotropic case at the same mag- 
netic field strength B, and <I>o is the quantum of flux. We 
look for this criterion to be significantly violated in one 
direction before the other. The factor of 1/2 allows the 
Lindemann parameter c to recover the usual definition in 
the isotropic case. 

For short wavelengths, the physics of an anisotropic su- 
perconductor can be mapped onto an isotropic supercon- 
ductor by a scaling procedure introduced by Blatter et al. 
|l3| Were this true at all wavelengths, there would be no 
reason to expect anisotropic melting to occur. However, 
scaling breaks down for the long wavelength bulk and 
tilt modes of the system, which are isotropic and have 
a vanishing energy cost. The result is that the spatial 
profile of the fluctuations of the vortices is less eccentric 
than the equilibrium lattice, suggesting an instability to 
partially melted (liquid crystalline) phases. 

Using the scaled momenta, q — (jk^/A, ky/^K, k^/A), 
the average fluctuations may be written as: 
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where the matrix C (q) and the elastic constants therein 
are functions of q, and the cutoff A — ^ °^ i ^"^^ 
by the vortex lattice spacing. The integrals are functions 
only of K, 7, and b = Tf^yy- We compute the integrals 
numerically to obtain the melting curves. 

We first compare to data on optimally doped YBCO, 
with i? ± c. It is believed that in this geometry, the in- 
trinsic pinning of the planes leads to a partially melted 
phase which is a smectic-C, in which vortices have melted 
along the planes. In the case of a lattice which is com- 
mensurate with the planes, the explicit symmetry break- 
ing of the planes adds a momentum-independent pinning 
term to the matrix element Cyy and tends to encour- 
age melting along the planes, into a smectic-C. In this 
sense, pinning competes with the aforementioned ten- 
dency of the anisotropy to encourage a smectic-A. We 
find for this material that the Lindemann criterion is vi- 
olated in one direction well before the other, as shown 
in Fig. ^, and the instability favors a smectic-C. For an 
incommensurate lattice, pinning is irrelevant, in which 
case the melting may proceed as in Fig. ^ as discussed 
below. For optimally doped YBCO with _B _L c, the 
lattice remains commensurate up to at least 120 Tesla, 
but in other types of superconductors it may be possi- 
ble to reach incommensurability at lower fields. jl5j Note 
that our approach is only capable of calculating the first 
melting curve, from solid to smectic. To calculate the 
melting curve for smectic to nematic, it is necessary to 
first derive elastic constants for the smectic phase, which 
is beyond the scope of the present paper. 

In Fig. ^, we plot the results in the absence of intrinsic 
pinning. The instability now favors a smectic-A, in which 
the long direction of each constituent "molecule" (in this 
case, the major axis of each cross sectional ellipse of a 
vortex) is oriented perpendicular to the melted layers, 
as in Fig. |[ Although this means that the vortices 
have melted first along the direction of harder shear, Cgg, 
this is the most common smectic geometry observed for 
oblong molecules. We show the schematic phase diagram 
in Fig. 0. The smectic region may be pinched off by first 
order transitions from lattice directly to nematic near 
Hc2{T = 0) and near T^. 

The problem of vortices in a 3D superconductor may 
be mapped to that of 2D bosons at zero temperature, 
by approximating the vortex interactions as "local" in 
the coordinate z, where the path of the vortex repre- 
sents a bosonic world line. The theory of anisotropic 2D 
melting in the presence of explicitly broken rotational 
symmetry predicts a region of quasi-long-range ordered 
(QLRO) smectic- A at finite temperature. [|l6[ ^ At zero 
temperature (to which the current case maps), a long 
range ordered (LRO) smectic-^ is possible. The inter- 
mediate smectic-A has also been seen in recent numerical 
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FIG. 3: Solid squares are vortex lattice melting data on op- 
timally doped YBCO with B _L c[Q. Open symbols are the 
results of the numerical integration of Eqn. |^, taking into 
account the pinning of the planes. Circles refer to melting in 
the "short" direction, and diamonds to melting in the "long" 
direction. We have taken the following parameters, appropri- 
ate for optimally doped YBCO with B L c : Tc = 92. 3A', 
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simulations of a two dimensional vortex system ||i8||. 

For the present case, the soft rotational modes usu- 
ally responsible for preventing translational LRO in the 
smectic are absent because the mass tensor introduces 
explicit rotational symmetry breaking. It costs energy to 
rotate the vortex smectic, and the system exhibits gradi- 
ent elasticity. It follows that a 3D smectic with explicitly 
broken rotational symmetry can have translational LRO. 
An interesting consequence of this is that the rigidity of 
the vortex smectic preserves superconductivity between 
the liquid-like layers, so that the transition from vortex 
lattice to vortex smectic is also a transition from 3D su- 
perconductivity to 2D superconductivity. 

Although we have presented results for a (homoge- 
neous) uniaxial superconductor, we also expect the re- 
sults to apply for a biaxial superconductor, with three 
different entries in the mass tensor. The cuprates cer- 
tainly exhibit anisotropy between the c direction and the 
planar directions, but they often also exhibit anisotropy 
within the a6-plane. In particular, our assumptions of 
mass anisotropy with no explicit translational symmetry 
breaking in the electronic degrees of freedom [pO[ are in 
principle satisfied for the geometry i3| |c in the cuprates in 
the presence of an electron nematic phase |2^ within the 
planes. Our assumptions may also be satisfied in stripe 
ordered phases, provided the mutual pinning is not too 
strong and thermal depinning of vortices from stripes oc- 
curs at a lower temperature than that at which the vortex 
lattice melts. 
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FIG. 4: Results of the numerical integration of Eqn. ^ in 
the absence of pinning. The circles refer to melting in the 
"short" direction, and the diamonds to melting in the "long" 
direction. We have taken the following parameters: 7^* — 
^ = 10, K = 1^ = 100, and H^2{T = 0) = lOOT. The 
figure is plotted for c — .2, as a function of & = ^ ^^^^ and 



The sniectic-j4 has clearest implications for experimen- 
tal probes that are capable of measuring the structure 
function of the vortex order, such as Bitter decoration 
(which is surface sensitive) or neutron scattering (which 
is a bulk probe). In melting from lattice to smectic, the 
diffraction pattern loses most Bragg peaks, retaining at 
most a line of Bragg peaks along the direction of modu- 
lated density, in this case a modulation of average mag- 
netic field density. More commonly, smectics exhibit the 
central Bragg peak along with a pair of peaks associated 
with the first harmonic of the density modulation. 

There are also distinctive implications of the double 
melting for ^SR. Muon spin rotation detects the dis- 
tribution of local magnetic fields. In the nematic phase, 
the time-averaged magnetic field density is uniform. In 
partialy freezing from the nematic to the smectic, fiSR 
would exhibit a new inhomogeneity in the magnetic field 
in the smectic state. Upon freezing further into the vor- 
tex lattice, the fiSR signal would reveal another transi- 
tion to further magnetic inhomogeneity. The changes in 
the fiSR signal are expected to coincide with the onset of 
highly anisotropic resistivity in going to the smectic, and 
with the onset of 3D superconductivity upon entering the 
lattice phase. 

Resistivity measurements are also sensitive to smectic 
order. When the Lorentz force is along a liquid-like di- 
rection, vortices move easily and the resistivity is large. 
When the Lorentz force is along the solid-like direction, 
the rigidity of the smectic resists vortex motion, and the 
resistivity vanishes (although the movement of defects 
may provide some small amount of dissipation). If the 
current is along the magnetic field direction, the Lorentz 
force and dissipation are negligible. The vortex smectic- 
A retains 2D superconductivity between the liquid-like 



layers of vortices, with the resistivity /5|| vanishing paral- 
lel to the layers, but p± ^ for currents applied perpen- 
dicular to the smectic layers. 

In conclusion, we have studied the problem of vor- 
tex lattice melting in anisotropic superconductors in the 
continuum limit. The introduction of anistropy in the 
mass tensor leads to elongation of vortex cross sections 
and interactions. We have demonstrated that interact- 
ing elongated vortices can form liquid crystalline phases. 
Using elasticity theory, with momentum-dependent elas- 
tic constants derived from Ginzburg-Landau theory, we 
have calculated the thermal fiuctuations of the vortex 
lattice. Comparing these results to an anisotropic Linde- 
mann criterion, we argue that there is an instability to 
an intermediate smectic phase. In the absence of intrin- 
sic pinning, we find an instability favoring a smectic-A, 
wherein the lattice has melted along the direction of the 
shorter lattice constant. 
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